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We study thermal instability in NbN superconducting stripline resonators. The system exhibits 
extreme nonlinearity near a bifurcation, which separates a monostable zone and an astable one. The 
lifetime of the metastable state, which is locally stable in the monostable zone, is measure near the 
bifurcation and the results are compared with a theory. Near bifurcation, where the lifetime becomes 
relatively short, the system exhibits strong amplification of a weak input modulation signal. We 
find that the frequency bandwidth of this amplification mechanism is limited by the rate of thermal 
relaxation. When the frequency of the input modulation signal becomes comparable or larger than 
this rate the response of the system exhibits sub-harmonics of various orders. 



Stochastic resonance (SR) is a phenomenon in which 
metastability in nonlinear systems is exploited to achieve 
amplification of weak signals [H, [3, 0] ■ SR has been ex- 
perimentally demonstrated in electrical, optical, super- 
conducting, neuronal and mechanical systems 0, H, H, 
0, i, S, Ea El, [II- Usually, SR is achieved by oper- 
ating the system in a region in which it has more than 
one locally stable (metastable) steady state and its re- 
sponse exhibits hysteresis. Under some appropriate con- 
ditions, a weak input signal, which modulates the transi- 
tion rates between these states, can lead to synchronized 
noise- induced transitions, allowing thus strong amplifica- 
tion. 

In the present paper we investigate SR and amplifica- 
tion in a superconducting (SC) NbN stripline resonator. 
Contrary to previous studies, we operate the system near 
a bifurcation between a monostable zone, in which the 
system has a single metastable state, and an astable zone, 
in which this state ceases to exist and the system lacks 
any steady states. In our previous studies we have in- 
vestigated several effects, e.g. strong amplification fl3| . 
noise squeezing jl^l , and response to optical illumination 



14l . Il5l | , which occur near this bifurcation, and limit cy- 



cle oscillations, which are observed in the astable zone 
T^ . 17 1 . In the present work we investigate experimen- 



tally and theoretically the response of the system to am- 
plitude modulated input signal, and find an unusual SR 
mechanism that has both properties of strong responsiv- 
ity and non-hysteretic behavior. The frequency band- 
width of this mechanism is found to be limited by the 
rate of thermal relaxation. We find that rather unique 
sub-harmonics of various orders are generated when the 
modulation frequency becomes comparable or larger than 
the relaxation rate. Moreover, we measure the lifetime 
of the metastable state in the monostable zone near the 
bifurcation and compare the results with a theory. 

Our experiments are performed using a novel device 
that integrates a narrow microbridge into a SC stripline 
electromagnetic resonator (see Fig. [I](A)). Design con- 
siderations, fabrication details as well as resonance modes 
calculation can be found elsewhere [l3| • The dynamics of 
our system can be captured by two coupled equations of 



motion, which are hereby briefly described (see Ref. [l3| 
for a detailed derivation) . Consider a resonator driven by 
a weakly coupled feed-line carrying an incident amplitude 
modulated coherent tone 5'" = 6q"(1 -|- a cos(wi„^))e~'"p', 
where 6™ is constant complex amplitude, Wp is the driv- 
ing angular frequency, a is the modulation depth, and 
ijjra ^ Wp is the modulation frequency. The mode am- 
plitude inside the resonator can be written as Be"*'^^*, 
where B (t) is a complex amplitude, which is assumed to 
vary slowly on a time scale of l/t^p- In this approxima- 
tion, the equation of motion of B reads [3] 



dB_ 
dt 
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where ujq is the angular resonance frequency and 7 (T) — 
7i + 72 (T), where 71 is the coupling coefficient between 
the resonator and the feed-line and 72 (T) is the temper- 
ature dependant damping rate of the mode, and T is the 
temperature of the microbridge. The term c™ represents 
an input Gaussian noise. The microbridge heat balance 
equation reads 



dT 9 



H{T- To) 



(2) 



where C is the thermal heat capacity, H is the heat trans- 
fer coefficient, and To = 4.2 K is the temperature of the 
coolant. 

Coupling between Eqs. (P) and ([2]) originates by the 
dependence of the damping rate 72 (T) of the driven 
mode on the resistance of the microbridge [l^ . which 
in turn depends on its temperature. We assume the sim- 
plest case, where this dependence is a step function that 
occurs at the critical temperature Tc ~ 10 K of the su- 
perconductor, namely 72 takes the value 72s for the SC 
T <Tc phase of the microbridge and 7211 for the normal- 
conducting (NC) T >Tc phase. 

Solutions of steady state response to a monochromatic 
excitation (no modulation a = 0) arc found by seeking 
stationary solutions to Eqs. H]) and ([2]) for the noiseless 
case c'" — 0. Due to the coupling the system may have, 
in general, up to two locally-stable steady-states, corre- 
sponding to the SC and NC phases of the microbridge. 
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FIG. 1: (A) Experimental setup. (B) System stability dia- 
gram. 
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The stability of each of these phases depends on both 

I ■ 1 2 

the power, Ppump oc \b™\ , and frequency Up parame- 
ters of the injected pump tone. Our system has four 
stabihty zones (Fig. [Hb)) [T^j. Two are mono-stable 
zones (MS(S) and MS(N)), where either the SC or the 
NC phases is locally stable, respectively. Another is a 
bistable zone (BiS), where both phases are locally stable 
13, [Hj. The third IS an astable zone (aS), where none of 
the phases are locally stable. Consequently, when the res- 
onator is biased to this zone, the microbridge oscillates 
between the two phases. The onset of this instability, 
namely the bifurcation threshold (BT), is defined as the 
boundary of the astable zone (see Fig. [Hb)). 

The experimental setup is depicted in Fig. [B^a). We 
inject an amplitude modulated pump tone into the res- 
onator and measure the reflected power in the frequency 
domain using a spectrum analyzer and in the time do- 
main using an oscilloscope. The parameters used for the 
numerical simulation were obtained as follows. The cou- 
pling coefficient 71 = 2 MHz and the damping rates 72s — 
2.2 MHz, 72n = 64 MHz were extracted from frequency 
response measurement 14 , 2^ , whereas the thermal heat 



capacity C = 54nJcm~ K~ and the heat transfer coef- 
ficient H = 12Wcm~^K~^ were calculated analytically 
according to Re fs. (2^. [23t. 

Our system exhibits an extremely strong amplifica- 
tion when tuned to the BT. Figure [2] shows both ex- 
perimental (Blue curves) and numerical (Red curves) 
results for the case where the system is driven by a 
modulated pump tone having the following parameters: 

0.0024, 
= 75 K. 



LOp = u)q — 2nx 4.363 GHz, uj^^ = 27r x 10 kHz, a - 
and by an effective noise temperature of T^g 



Panel (A) plots the signal gain Ggig, defined as the ratio 
between the reflected power at frequency Wp-t-Wm and the 
sum of the injected powers at frequencies Up ± w,n, as a 
function of the mean injected pump power (Ppump). The 
system exhibits large gain of approximately 20dB around 
the BT. The experimental results exhibits excess gain be- 
low BT relative to the numerical results. This can be 



explained by additional nonlinear mechanisms 2J] that 
may induce small amplification, and are not theoretically 
included in our piecewise linear model. 

Figure[2](B), shows time and frequency domain results 
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FIG. 2: (A) Experimental (dotted-blue) and numerical 
(crossed-red) results of the signal amplification Gsig as a func- 
tion of the mean injected pump power < Ppump >. (B) Ex- 
perimental (subplots (xV) and (2:2)) and numerical (subplots 
(a;3) and (a;4)) results of the reflected power Prcfl as a function 
of time (subplots (a;l) and (a;3)) and scanned frequency /sa 
(subplots (a;2) and (a;4)), centered on the resonance frequency 
/o — 4.363GHz (A/sA ~ /sa — /o), where x denotes a, &, and 
c, corresponding to the marked points in panel (A). The 
dashed-green curve represents the modulation signal. The 
time domain measurements are normalized by their maximum 
peak to peak value. 



of the reflected power, for three pairs of input power val- 
ues, corresponding to the marked points (a — c) in panel 
(A) . In addition, the time domain measurements contain 
a green-dashed curve showing the modulating signal. The 
results shown in subplots (al — a4) were obtained while 
biasing the system below the BT, namely, (Ppump) was 
set below the power threshold, P^- In general, the spikes 
in the time domain plots of Fig. [2] (B) indicate events 



in which the temperature T temporarily exceeds 17 1 . 
Below threshold, the average time between such events, 
which are induced by input noise, is the lifetime F"-'^ of 
the metastable state of the resonator. As we will show 
in the last part of this paper, F strongly depends on the 
pump power near BT, thus power modulation results in 
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a modulation of the rate of spikes, as can be seen both 
in the experimental and simulation results. 

Subplots (bl — b4) of Fig. [DJB) show experiments 
in which (Ppump) — and thus, the modulation itself 
drives the resonator in and out the astable zone. As a 
result, during approximately half of the modulation pe- 
riod nearly regular spikes in reflected power are observed, 
whereas during the other half only few noise-induced 
spikes are triggered. This behavior leads to a very strong 
gain as well as to the creation of higher order frequency 
components (subplots (b2,b4)). Figure [2jB), Subplots 
(cl — c4), show experiments in which (Ppump) > ^c, and 
thus the regular spikes occur throughout the modula- 
tion period. The rate of the spikes is strongly correlated 
to the injected power [l^, and it is higher for stronger 
pump powers. Therefore, as the injected pump power is 
modulated, so is that rate. This behavior also creates a 
rather strong amplification, though weaker than the one 
achieved in the previous case. 

The amplification mechanism in our system is unique 
in several aspects. First it is extremely strong. To 
emphasize the strength of the amplification we note 
that, usually, no amplification greater than unity (0 dB) 
is achieved in such measurements with SC resonators 



26l |. unless the resonator is driven near BT [27|. In 
addition, it does not exhibit a hysteretic behavior. 

Each spike in subplots (al — a2) of Fig. lasts ap- 
proximately 1 ^s, after which the device is ready to detect 
a new event. This recovery time determines the detection 
bandwidth. A measurement of the dependence of the am- 
plification mechanism on the modulation frequency 
has reveled a mechanism in which sub-harmonics of the 
modulation frequency are generated by the device. The 
generation occurs when the modulation period is compa- 
rable to the recovery time of the system. The results are 
shown in Fig. [3] which shows both experimental (Blue 
curves) and numerical results (Red curves) for the case 
of Wp = 27r X 4.363 GHz, a = 0.017, T^s = 75 K, and 
cj,„ 27r X 1.2 MHz for panel (A) and uj^ 27r x 2 MHz 
for panel (B). Panel (A), shows the reflected power, ob- 
tained for three gradually increased pump power values, 
and corresponding to sub-harmonics generation (SHG) 
of the second, third, and forth orders. SHG of the third 
order, for example, are generated by a quasi-periodic re- 
sponse of the system (subplots (al,a3)). Each quasi- 
period lasts three modulation cycles, where only during 
the first two a spike occurs, namely a spike is absent 
once every three modulation cycles. This behavior origi- 
nates from the mismatch between the modulation period 
and the recovery time of a spike, which induces a phase 
difference, that is monotonically accumulated, between 
the two. Once every n = 3 modulation cycles, in this 
case, the system fails to achieve critical conditions near 
the time where the peak in the modulation occurs, and 
therefore a spike is not triggered. Similar behavior is also 
shown in subplots (al,a3) and (cl,c3), where the quasi- 
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FIG. 3: Sub-harmonics generation. 

period lasts two and four modulation cycles respectively. 

Another mechanism for SHG is observed when the 
modulation frequency is increased. Fig. [3l Panel (B), 
shows measurement results for = 2tt x 2 MHz, which 
demonstrate SHG of order n — 3. Unlike the previous 
case, this SHG is characterized by a single spike that 
occurs once every three modulation cycles. 

We further study our system by measuring the 
fluctuation-induced escape rate F of the metastable state 
in the MS(S) zone. In Ref. we have found theoreti- 
cally that 



F = Fo exp( 



71 AP2 

' ^ pump 



(3) 



pump 



where Fq = y/HjJC /2tt, and the power difference is 
given by APpump = Pc — Ppump- Note that the unusual 
scaling law in the present case log(F/Fo) oc APp^j^^p, 
which differs from the commonly obtained scaling low 
of log(F/Fo) (X APpump m, Is^l, is a signature of the 
piecewise linear dynamics of our system. 

The escape rate was experimentally measured for sev- 
eral levels of Teff, which are given in the first row of 
table m The noise was generated by an external white 
noise source, and combined with the amplitude modu- 




FIG. 4: Escape rate of metastable states for several levels of 
Tcff, summarized in table HI The graphs are plotted in pairs, 
where the solid curves show the experimental data and the 
dashed curves show the corresponding theoretical fit. 



lated pump tone. The modulation frequency was set 
to 500 Hz, which is more than three orders of magni- 
tude lower than the relaxation rate of the system, and 
therefore to a good approximation the system follows this 
modulation adiabatically [29|]. 

The results are shown in Fig. 31 which plots the escape 
rate in logarithmic scale as a function of APp^^^p. Six 
pairs of solid and dashed curves are shown, correspond- 
ing to the six different levels of injected noise intensities. 
The solid curves were extracted from time domain mea- 
surements of the reflected power. The dashed curves were 
obtained by numerically fitting the experimental data to 
Eq. ^ and show good quantitative agreement between 
the experimental results and Eq. The fitting param- 
eters included the pre-factor Fq — 0.86 MHz that was 
determined by a separate fitting process, and Pc (see 
table m that slightly decreases with the thermal noise. 
This behavior can be explained by local heating of the 
microbridge, induced by the noise that is injected into 
the resonator through additional resonance modes. Note 
that Toff was extracted from a direct measurement of the 
injected noise intensity (see table |T|. Note also that the 
system recovery time at the threshold imposes a limit on 
the measured escape rate. Thus the escape rate close to 
the threshold might be higher than measured. 

In summary, a novel mechanism of SR with a single 
metastable state has been demonstrated. Near BT the 
system exhibits rich dynamical effects including bifurca- 
tion amplification and SHG. In spite of its simplicity, our 
theoretical model successfully accounts for most of the 
experimental results. 
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TABLE I: Escape rate parameters 





1 


2 


3 


4 


5 
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Note 


Tcff [10" K] 


0.52 


1 


1.36 


1.64 


2.3 


3.76 


Measured 




125.1 


121.6 


120.5 


120 


119.2 


117.6 


Fitted 
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